Emergent motion of condensates in mass-transport models 
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We examine the effect of spatial correlations on the phenomenon of real-space condensation in 
driven mass-transport systems. We suggest that in a broad class of models with a spatially correlated 
steady state, the condensate drifts with a non-vanishing velocity. We present a robust mechanism 
leading to this condensate drift. This is done within the framework of a generalized zero-range 
process (ZRP) in which, unlike the usual ZRP, the steady state is not a product measure. The 
validity of the mechanism in other mass-transport models is discussed. 
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Nonequilibrium condensation, whereby a macroscopic 
fraction of microscopic constituents of a system accumu- 
lates in a local region, is a common feature of many mass- 
transport systems. Examples include shaken granular 
gasses [T], vehicular traffic PHll, the macroeconomics of 
wealth distribution [H [6] , and others [3 18] . Mechanisms 
which can lead to the formation of condensates have been 
studied extensively in recent years, mainly by analyzing 
prototypical toy models. A primary role in these studies 
was played by the zero-range process (ZRP), an exactly- 
solvable model in which particles hop between sites with 
rates which depend only on the number of particles in 
the departure site Extensions and variations of 

the ZRP have been used to study the emergence of multi- 
ple condensates [T^] , first order condensation transitions 
[131 [14] and the effect of interactions [T^ and disorder [16] 
on condensation. Moreover, one-dimensional phase sepa- 
ration transitions in exclusion processes and other driven 
diffusive systems can quite generally be understood by a 
mapping on ZRPs [I7j . 

The dynamics of condensates is less well explored. In 
the ZRP, where condensation takes place when a macro- 
scopic fraction of particles occupies a single site, the re- 
sulting condensate does not drift in the thermodynamic 
limit [131 [THHin] . It is shown below that this is related 
to the fact that the steady state of the ZRP is a product 
measure. In some real-world systems, however, conden- 
sates are in continual motion. For example, traffic jams, 
which can be viewed as condensates, are known to prop- 
agate along congested roads [TTJ [HI [35] . Recently, two 
variants of the ZRP were also found to relax to a time 
dependent state in which the condensate performs a drift 
motion: one is a ZRP with non-Markovian hopping rates 
[23l[23], and the other is a model with "explosive conden- 
sation" [5S]. To date, there is no systematic understand- 
ing of the mechanism by which a macroscopic condensate 
motion emerges from the underlying nonequilibrium mi- 
croscopic dynamics. 

In this Letter, we study how spatial correlations in 
the steady states may lead the condensate to drift with 
a non-vanishing velocity. We do so by introducing a 



generalization of the ZRP whose steady state does not 
factorize. Within its framework, we identify the mech- 
anism which generates the drift. The analysis is based 
on numeric simulations and on a mean-field approxima- 
tion which captures the essential effect of correlations in 
the condensed phase, and thus elucidates the different 
observed modes of condensate motion. The drift mech- 
anism which we identify is robust and therefore it is ex- 
pected to be valid in a broad class of spatially-correlated 
mass-transport systems. 

We focus on a class of stochastic one-dimensional mod- 
els defined on a ring with L sites. At any given time, each 
site i is occupied by particles with — N, rij > 0. 

The model evolves by a totally asymmetric hopping pro- 
cess whereby particles hop from site i to i + 1 with a rate 
which depends on the occupation numbers rii and n^-i. 
This is a generalization of the usual ZRP in which the 
rate depends only on n^. More specifically, we choose the 
hopping rates to be of the form 

■w{ni-l)u(ni) 

n,_i,ni,ni+i > ni_i,n, - l,rii+i -I- 1, (1) 

with rates 

The particular form of u{n) is motivated by the fact that 
in the usual ZRP, which corresponds to a = 1, this choice 
with b > 2 leads to a condensation transition [26]. The 
rate w with a ^ 1 represents an interaction between 
nearest-neighbor sites. According the dynamical rules 
([T])-([2]), at every short time interval dt, each site i whose 
occupation Ui > 1 may eject a particle with a probability 
u{ni)dt, as long as the preceding site (i — 1) is occupied. 
If the preceding site is empty, this probability changes to 
au{ni)dt. The model has three parameters: 5, a, and the 
density p = N/L which is conserved by the dynamics. 

In the usual ZRP (the case of a = 1), the stationary 
distribution is known to factorize into a product of single 
site terms, and so can be exactly calculated [H [27l l28] . 
This factorization property renders the ZRP quite spe- 
cial, as slight variations of the ZRP dynamics result in 
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FIG. 1. The location (top) and occupation (bottom) of the most occupied (•) and 2"^^ most occupied (o) sites for several values 
of a. (a) a = 1, i.e., the usual ZRP. The condensate is stable for long times, and relocates to random distant sites, (b) a = 1.5. 
The condensate advances through a slinky motion from one site to the next, (c) a = 1.05. The condensate remains on each site 
for a long time before "spilling" to the next. The definitions of Tbarrier and Tspiii are indicated, (d) a = 0.5. The condensate 
skips every other site. In all cases 6 = 3, p = 10 and L = 1000, except (a) where L = 400. Note the different time scales. 



non-factorizable models. To probe the effect of spatial 
correlations on the condensate, we choose for simplicity 
w to be of the form ([2]), which leads to a spatially corre- 
lated steady-state when a ^ 1. The drift motion which 
we describe below is also found for other forms of w(n), 
such as cases where w{n) ^ 1 for finitely many values of 
n, and for other choices of u{n) which give rise to con- 
densation 29J. 

We have carried out Monte Carlo simulations of the 
dynamics ([T])-([2| for several values of a in a system of 
size L = 1000 and density p = 10. After the system has 
relaxed to its steady state, the dynamics of the conden- 
sate was examined by tracking the position of the most 
occupied sites over time. The results are presented in 
Fig. [I] and in videos in the supplemental material [30] . 
In the usual asymmetric ZRP (Fig.jTJi), it is known that 
the condensate is static up to timescales of order and 
then it relocates to a random site due to fluctuations 
[HI [T51 - BD] . There is a striking qualitative difference in 
the dynamics of model when a^l (Fig. [l]D-{lJi), 

where the condensate is clearly seen to drift along the 
lattice. The condensate is seen to move from one site to 
the next when a > 1 (Fig. [T|d), or to skip every other site 
when a < 1 (Fig. [iji). In both cases, when a is not too 
close to 1 the motion is "slinky" -like, with the conden- 
sate "spilling" from an old site to a new one immediately 
after the previous spilling has completed. The drift be- 
comes somewhat less regular when a is close to 1. In 
this regime, the slinky motion is interrupted by periods 
of time when the condensate occupies a single site, before 
the spilling process is initiated (Fig. This, however, 
is argued below to be a crossover mode, and the interval 
in a in which it is observed shrinks in the large L limit. 

To understand these results we propose a mean-field 
analysis of the model in which the occupations of all 
sites are considered independent, but might not be identi- 
cally distributed. Within this approximation, the current 
which arrives into site i from site i — 1 is a Poisson pro- 
cess whose rate we denote by J;. The probability Pi{ni) 



to find Ui particles in site i thus evolves according to 
dP-in) 

= P^{n^ - ^)J^ + P^{n^ + l){w^)u{n, + 1) 

^P,{ni){j, + {wi)u{n,)), (3) 

where 

oo 

{w,) = J2 P^^l{n)w{n) = !-!-(«- l)P.-i(O) (4) 

n=0 

encodes the mean effect of site i 1 on the hopping rate 
out of site i. Equation ([3| is valid also when Ui — with 
the definition Pi(— 1) = 0. Equations ([s]) and (|4| are 
to be solved with the self consistency condition J^+i ~ 

At low density, the system is in a subcritical, disor- 
dered phase (this will be shown below). In this homo- 
geneous phase, Pi{n) = P{n) and Ji = J for all sites i. 
At higher densities, however, condensation takes place, 
where the translational symmetry is spontaneously bro- 
ken and both Pi (n) and Ji depend on the distance of site i 
from the condensate. This dependence of P and J on i is 
a result of the correlations which exist in the steady state 
of the model, and it provides the mechanism for the con- 
densate drift: a nonhomogeneous Ji implies that in some 
sites the outflowing current is smaller than the incoming 
current, leading these sites to accumulate particles while 
other sites are similarly being depleted of particles. We 
shall now demonstrate that this occurs in our model. 

In the homogeneous (subcritical and critical) phases, 
the model eventually reaches a steady state. In the non- 
homogeneous supercritical phase, however, the conden- 
sate location keeps moving with time. The analysis of 
this time-dependent phase is based on one key observa- 
tion: the timescale of the microscopic dynamics, which 
for the rates ^ is of order 1, is much faster than that of 
the condensate motion. As shown below, the timescale of 
the spilling process is of order L, validating this observa- 
tion in the thermodynamic limit. Due to this timescale 
separation, while the condensate (i.e. the most occupied 



3 



site) , is static all other sites reach a quasi-stationary dis- 
tribution. 

In both phases, by equating the LHS of Eq. (|3| to zero 
the (quasi-)stationary distribution is found to be 



with Zi = Ji/{wi). (5) 



k<r 



Here, Zi plays the role of a "fugacity" of site i. For rates 
of the form ([2|, the normalization of Pi{n) yields 



Pm = [2^1(1,1; 6 +l;z,)]" 



(6) 



where 2^1 is a hypergeometric function (note that P,;(0) 
depends on the exact form of u{n) and not just on its 
large n asymptotics) . The occupation probability is 
asymptotically given by Pi{n) ~ n~^zf . 

We first examine the solution ([s]) in the subcritical and 
critical phases, and show that the model undergos a con- 
densation transition. Since the system is homogenous in 
these phases, the subscript i may be dropped from equa- 
tions (|4])-([6|. The fugacity can now be determined in 
terms of the density by inverting the relation 

, . p. N 2^1 (2, 2; 6 + 2; z) 

^ (l + 6)2Fi(l,l;&-|-l;z) 

where the RHS is obtained by substituting Eq. ([s]) in the 
sum. Similarly, Eq. Q for {w) reads in the homogeneous 
phases (w) = 1 + (a - l)[2Fi(l, 1;6 + 1;2:)]"^ 

The density Q is an increasing function of z that at- 
tains its maximum at z = 1, which is its radius of con- 
vergence about the origin. A finite density at z = 1 
indicates a condensation phase transition, which is math- 
ematically similar to Bose-Einstein condensation By 
substituting z = 1 in ([?]) it is seen that condensation 
takes place when h > 2, in which case the critical den- 
sity is pc — 1/(6 — 2), the same value as that of the 
usual ZRP. The critical current is similarly found to be 
Jc = {w)z^i = 1 -I- (a — i)b/{b — 1). As long as p < Pc, 
the system remains in a homogeneous subcritical phase. 
When p is increased, the current J increases until p and 
J reach their critical values and all sites of the system 
are in a homogenous critical phase. When the density 
is further increased, condensation sets in, breaking the 
translational invariance of the system. 

Let us now discuss the nonhomogeneous supercritical 
phase and the mechanism of the condensate motion. We 
focus on the case of a > 1. In this case, the condensed 
phase is composed of a condensate, which at any given 
time consists of two macroscopically occupied consecu- 
tive sites (say 1 and 2), while the rest of the sites are mi- 
croscopically occupied. We show that J2 > 1 and Ji = I 
for z ^ 2. This results in an increase of the occupa- 
tion of site 2 at the expense of site 1 over a macroscopic 
0{L) time scale, while the rest of the sites are in a quasi- 
stationary state. Therefore, the condensate drifts with a 
velocity of order L~^. 



The analysis begins at site 1, whose occupation we as- 
sume is ni = 0{L) ^ 1, and thus it emits a mean current 
J2 — (wi)(l + {b/ni)) ~ {wi). At the moment, (wi) is 
unknown. It is determined self-consistently at the end of 
the calculation. Since Pl{0) 7^ 0, as is established below, 
it is seen that J2 > 1 (since a > 1). We now proceed 
to examine the second site. As long as site 1 accommo- 
dates the condensate it is never empty, i.e., Pi(0) — 0. 
It follows from Q that {W2} — 1- The fugacity of the 
second site is then Z2 = ^2/(^2) — >^2 > 1, and there- 
fore its occupation distribution ^ cannot be normalized. 
This means that as long as site 1 is highly occupied, site 
2 tends to accumulate particles, implying that its occu- 
pation too becomes macroscopic (of order L) for a long 
period of time [3TJ [32] . We call such a site with fugacity 
Zi > 1 supercritical. 

The analysis now continues site by site in a simi- 
lar fashion. For each site i, (wi) is calculated using 
Q from the known value of Pj_i(0). The fugacity 
of site i ([5]) is then calculated from (wi) and the in- 
coming current into the site Ji. Once the fugacity is 
known, Pi{0) and Ji+i are determined from (|6| and from 
Ji_i_i = X^n ^^'^ process is repeated 
in the next site. Performing this analysis reveals that 
site 3 is critical (i.e., Z3 = 1) and sites i — 4, 
are subcritical (z^ < 1), with Jj = 1 and (wi+i) — 
l + {a-l)/2Fi{l,l;b+l]l/{w,)) for aU i > 3. This re- 
cursion relation defines a sequence (wi) which converges 
(exponentially rapidly) to a unique fixed point w* (a, b) 
which is the solution of the equation 



(8) 



2Fi(l,l;5+l;l/w;*)' 



and thus satisfies w*{a,b) > 1 for all a > 1. When 
L ^ 1, the periodic boundary conditions imply that 
{wi) ~ w* > 1, and thus Eq. ^ confirms that Pl{0) > 0. 
This closes the loop self-consistently and completes the 
calculation of the quasi-stationary distribution for the 
non-homogenous phase. 

A natural order parameter for the condensation tran- 
sition is the bulk density of the "background fluid" pbg, 
which can be defined as the mean density of all but 
the two most occupied sites (since the condensate is 
typically carried by two sites). Below the transition, 
Pbg = P, which approaches pc — l/{b — 2) as the tran- 
sition is approached from below. Above the transition, 
all sites outside of a finite boundary layer around the 
condensate are subcritical with a mean occupation of 
Pbg — p{z = ^/w*) < p{l) = Pc since the function p(z), 
Eq. ([7]), increases monotonically with z. Therefore, the 
condensation transition is found to be a discontinuous 
(first order) one. This is in contrast to the usual ZRP 
with rates ^ and a = 1 where the transition is continu- 
ous. A similar discontinuity exists in the current, which 
jumps from Jc > 1 just below pc to J = 1 just above it. 
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We now discuss the emergent dynamics of the conden- 
sate and identify two distinct modes of motion: a reg- 
ular shnky motion, and an irregular motion through a 
barrier. The motion of the condensate from one site to 
the next consists of two stages: a "spilling" stage during 
which it is supported on two sites, and a period before 
this spilling is initiated, when the condensate is carried 
by a single site. We first consider the spilling process. 
According to the calculation above, the number of par- 
ticles that accumulate in the second condensate site per 
unit time is on average J2 — J3 = w* — 1. As there are 
A^cond — {p ~ Pbg)L particles in the condensate, the to- 
tal spilling time Tgpin scales, to leading order, linearly 
with the system size: Tspni = {p — pbg)L/{w* — 1). This 
justifies the assumption of timescale separation which un- 
derlies the existence of a quasi-stationary state. In the 
limit of a — >■ 1, the spilling time diverges. 

Once a spilling is complete, there is a moment that the 
condensate is located solely on a single site. We now re- 
label this site as site 1. At this moment, the occupation 
of the following site is ri2 ~ Pbg — 0(1). The rate at 
which particles leave the second site is, at this stage, ap- 
proximately J3 « 1 -I- 5//5bGj which should be compared 
with the rate of incoming particles, J2 — w* . According 
to Eq. ([s]) and the definition of pbg, the two rates are 
equal when a = a* which is the solution of the equation 
a* = l + b2Fi{l,l;b+l;l/w*{b,a*))/pBG{b,a*). The 
mode of condensate motion now depends on whether a 
is larger or smaller than a* . (i) When a > a* , the initial 
current into site 2 is larger than the mean current out 
of this site, and a spilling of the condensate is initiated 
immediately. In this case, the condensate drifts continu- 
ously in a slinky motion as in Fig. [IJj. (ii) On the other 
hand, J2 < J3 when 1 < a < a*, and thus particles 
do not immediately accumulate on site 2. The incoming 
current into the site surpasses the outgoing current and 
spilling sets in only after fluctuations bring the occupa- 
tion of the second site to a value n* (a, b) which is defined 
by ui* = 1 + b/n* . The ensuing motion of the condensate 
is more irregular, with a stable condensate which occa- 
sionally spills to the next site as in Fig. [T];. Note that the 
mean time Tbarrior it takes before a fluctuation brings 7^2 
over the barrier n* does not scale with the system size. 
Thus, in a large enough system the condensate regularly 
drifts and is typically supported by two neighboring sites 
for any value of a > 1. 

The numerical simulations support the qualitative pic- 
ture which emerges from the MF analysis presented 
above, even if the two differ quantitatively. The exis- 
tence of two modes of motion, slinky-like and irregular, 
and the crossover between them as a is increased con- 
form with numerical findings (Fig. [T]). In particular, the 
spilling mechanism between the two condensate sites in 
which the accumulation of particles is linear in time is 
verified (Fig.jlja). Furthermore, the first order nature of 
the transition, as manifested by the Pbg(p) curve, and the 




FIG. 2. A first order phase transition is seen in the back- 
ground density pbg as a function of the density p. Numerical 
results for several system sizes are plotted, along with an ex- 
trapolation to L = 00 (thick black line). Here q = 1.5 and 
6 = 3. The inset shows that the occupation of a site far away 
from the condensate has a subcritical distribution (i.e., with 
an exponential tail) when a > 1. This difi'ers from the usual 
ZRP where the background fluid is known to be critical. Re- 
sults are for L — 1000, p — 10, and 6 = 3. 

subcritical nature of the background fluid are presented 
in Fig. [2} 

The mechanism for the condensate drift found in this 
model can be summarized as follows: the spontaneous 
breaking of translation invariance by the formation of 
the condensate may induce an accumulation of parti- 
cles in a nearby site. This accumulation results in a 
continually drifting condensate, since whenever a con- 
densate is established on a new site, another one be- 
gins to form further ahead. This mechanism holds in 
a much more general setting, including when other forms 
w{n), partially asymmetric hopping and higher dimen- 
sional lattices are considered, and more widely in other 
non-factorized mass-transport models [22 ■ Note that the 
drift discussed here, in which the two most occupied sites 
are typically nearest neighbors, cannot occur in models 
with a factorized steady state, the latter being symmet- 
ric under site permutations. In this respect, our mech- 
anism differs from that studied recently in where 
unbounded hopping rates generate a drift (with infinite 
velocity) in a model whose steady state factorizes. 

An important point to note is that in general, the new 
condensate site does not have to be a neighbor of the old 
one. For instance, in our model ([T])-([2]) with a < 1, a sim- 
ilar analysis shows that the condensate skips every other 
site, as observed in Fig. [TJi [29j. In this case, the super- 
critical site is site 3, rather than 2 (when the condensate 
is located on site 1). In principle, it may happen that 
there is more than one supercritical site, possibly lead- 
ing to more complicated condensate drifts. It may also 
happen that no other site is supercritical, in which case 
the condensate would not drift. A precise and general 
classification of the conditions under which a condensate 
drift occurs remains an interesting open problem. How- 
ever, in many specific models, a study of condensation 
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and the condensate motion can be carried out following 
the mean-field procedure outlined in this Letter. For in- 
stance, a recently proposed accelerated exclusion process 
(AEP) [33] can be analyzed in a similar fashion, yielding 
the phase diagram of the model and revealing that the 
AEP condensate drifts in the steady-state [33]. 
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